arXiv: 1509.05425vl [astro-ph.GA] 17 Sep 2015 


Mon. Not. R. Astron. Soc. 000,[7|fT2l(2015) Printed 21 September 2015 (MN RTbX style file v2.2) 


BAT AGN spectroscopic survey II: X-ray emission and high 
ionization optical emission lines 

Simon Berney^, Michael Koss^^t, Benny Trakhtenbrot^l, Claudio Ricci^, Isabella 
Lamperti\ Kevin Schawinski\ Mislav Balokovic^, D. Michael Crenshaw"^, Travis 
Fischer^, Neil Gehrels^, Fiona Harrison^, Yasuhiro Hashimoto^, Kohei Ichikawa^, 
Richard Mushotzky^, Kyuseok Oh\ Daniel Stern^, Ezequiel Treister^^, Yoshihiro 
Ueda^, Sylvain Veilleux^, and Lisa Winter^^ 

^ Institute for Astronomy, Department of Physics, ETH Zurich, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland 
^Institute de Astrofi'sica, Facultad de Fi'sica, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 
^Cahill Center for Astronomy and Astrophysics, California Institute of Technology, Pasadena, CA 91125, USA 

‘^Department of Physics and Astronomy, Georgia State University, Astronomy Offices, One Park Place South SE, Suite 700, Atlanta, GA 30303, USA 
^NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

^Department of Earth Sciences, National Taiwan Normal University, No. 88, Sec. 4, Tingzhou Rd., Wenshan District, Taipei 11677, Taiwan R.O.C 
''Department of Astronomy, Kyoto University, Kyoto 606-8502, Japan 
^Astronomy Department, University of Maryland, College Park, MD, USA 

Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, MS 169-21, Pasadena, CA 9II09, USA 
Departamento de Astronomia, Universidad de Concepcion, Concepcion, Chile 
''Atmospheric and Environmental Research, 131 Hartwell Ave No. 4, Lexington, MA 02421, USA 
^Ambizione fellow 
^Zwicky fellow 


ABSTRACT 

We investigate the relationship between X-ray and optical line emission in 340 nearby 
(z ~ 0.04) AGN selected above 10 keV using Swift BAT. We find a weak correlation be¬ 
tween the extinction corrected [Olll] and hard X-ray luminosity (Tjoni] ^ 14 - 195 ) with a 
large scatter (Rpear = 0.64, O’ = 0.62 dex) and a similarly large scatter with the intrinsic 2—10 
keV to [O III] luminosities (Rpear = 0.63, O = 0.63 dex). Correlations of the hard X-ray fluxes 
with the fluxes of high-ionization narrow lines ([O III], He II, [Ne III] and [Ne v]) are not sig¬ 
nificantly better than with the low ionization lines (Ha, [S ll]). Factors like obscuration or 
physical slit size are not found to be a significant part of the large scatter. In contrast, the opti¬ 
cal emission lines show much better correlations with each other (ff = 0.3 dex) than with the 
X-ray flux. The inherent large scatter questions the common usage of narrow emission lines 
as AGN bolometric luminosity indicators and suggests that other issues such as geometrical 
differences in the scattering of the ionized gas or long term AGN variability are important. 
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1 INTRODUCTION 

Based on the unification model of active galactic nuclei (AGN; 
|Urry et ^ |1991[ |AntonuccT| |1993[ |Urry & Padovani] |1995| ), 
narrow optical emission lines are emitted from the narrow-line 
region (NLR). The NLR is ionized by the Extreme UV-soft 
X-ray emission from the AGN’s central engine. Afterwards, the 
illuminated interstellar gas in the NLR radiatively cools down by 
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emitting the narrow optical lines. Unlike this optical emission, 
X-rays are emitted from much closer to the supermassive black 
hole (SMBH). The UV photons radiated by the accretion disk are 
inverse Compton scattered inside the corona and can reach energies 
up to a few hundreds of keV (e.g., |Netzer|[20T^ . Thus, if AGN 
are accreting at a constant rate and constantly ionizing the NLR, 
we would naively expect a tight correlation between the X-ray 
emission from the corona and the optical emission lines from the 
NLR. This is why both high ionization optical emission lines and 
X-rays are thought to be reliable tracers of the AGN bolometric lu¬ 
minosity (L/,o/) and of the black hole growth ({Heckman et al.|2005^. 
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The bolometric luminosity {Lhoi) of unobscured AGN esti¬ 
mated from the spectral energy distribution (SED) fitting technique 
using optical, UV, and X-ray data shows a tight correlation with 
hard X-ray data alone (e.g., 14—195 keV with a scatter of 0.3 
dex, |Vasudevan et al.|2009^ though X-ray surveys may not detect 
heavily obscured Compton-thick AGN (Ah > 10^^ cm^^). Strong 
observational correlation has also been found between UV emis¬ 
sion and X-ray emission ( [Steffen et al.|200^ [Young et al.|2010| (. 
The use of [O III] A 5007 to infer AGN bolometric luminosities has 
long history (e.g., Heckman et al.[[2004l [LaMassa et al.[[200^ 


[Kauffmann & Heckman 20091. It has been widely used to estimate 
the accretion rates of AGN using bolometric corrections based on 
the X-ray emission, the 5100 A monochromatic continuum or the 
mid-IR emission (e.g., [Lamastra et al.[|2009[ . However, several 
studies have noted that the strength of [Olll] A5007 is modulated 
by factors including covering factor of the narrow-line region 
gas and its de nsity (e.g., [Net zer & Laor[ [1993) , and ionization 
parameter (e.g., [Baskin & Laor|2005^^ 


Previous studies have compared X-ray emission with optical 
emission lines strengths in hard X-ray selected samples. |Heckman| 
[et al.[ j2005[ used a sample of 47 local AGN (z <0.2) from the 
RXTE slew survey to analyze the correlation between 3-20 keV 
flux and the [O III] strength ^Revnivtsev et al.[[2004[ [Sazonov ^ 
[Revnivtsev[[2004|l. Other works have also made use of the early 
results from the Swift BAT survey (9-month catalog, [Tueller et al. 


20081 to analyze optical and X-ray emission from AGN. 


et al. 


Melendez 


j2008^ studied the correlation between the [O III] A5007 line 


strength and two X-ray bands (2—10 keV and 14—195 keV) for 40 
low redshift (z<0.08) BAT AGN from the 9-month catalog. [Wint er] 
[et al.[ ( |2010[ made a similar analysis with a slightly larger sample 
(64 BAT AGN). Finally, an analysis of the full 9-month catalog 
has recently been completed (sslOO BAT AGN, Ueda et al. 2015, 
submitted). [Netzer et al.[ ( |2006[ l studied a sample of 52 obscured 
AGN at high redshift (z ~1) suggesting that [Olll] emission is 
weaker relative to 2—10 keV emission in high Lx sources. 


2 SAMPLE 

The optical data are taken from the BAT AGN Spectroscopic Sur¬ 
vey (BASS) Data Release 1 (Koss et al., in prep.), which gathered 
559 optical spectra from targeted spectroscopic campaigns on BAT 
sources and public optical surveys (SDSS and 6DF, [Alam et al.[ 
|20T5l |7on^ et al.|2009[ (. This sample is the largest and the most 
complete catalog of AGN optical spectra selected from the Swift 
hard X-ray survey. It contains 67.6% of the total AGN detections 
in the Swift BAT 70-month catalog and has a mean redshift of 
z = 0.10. 

From this parent sample, we first remove the spectra coming 
from the 6dF survey because they lack flux calibration. We remove 
all blazars listed in the Roma blazar catalog (BZCAT) v5.0 
( [Massaro et al.|200^ since the X-ray emission in such sources is 
affected by differential beaming. We also apply two redshift cuts. 
The first one, z > 0.01, avoids very nearby star-forming galaxies 
(such as M82) and ensures that a significantly large portion of the 
ionized region within the galaxy is covered (i.e., >0.2 kpc). We 
also limit to z < 0.4 to avoid redder emission lines such as Ha 
from being shifted outside the observed wavelength interval. Addi¬ 
tionally, we remove ten spectra due to calibration problems around 
the [OIII] line. Finally, we exclude SDSS J155334.73-I-261441.4, 
a gravitationally lensed quasar ( [Ofek et al.[2007[ (. The final sample 
totals 340 optical spectra with (z) = 0.05. 

We define an emission line detection when we reach a SNR>3 
of the line with respect to the noise of the adjacent continuum. Out 
of these 340 objects, we have 327 detections of [O III], 301 of Ha, 
281 of [Nil], 278 of [S II], 265 of Hj3, 206 of [O I], 144 of [Ne III], 
106 of Hell and 51 of [Ne V]. The intrinsic 2—10 keV flux values 
and the column density measurements provided for 335/340 AGN 
are based on a homogeneous model fitting using the best available 
X-ray data with simultaneous modeling of the 0.2—10 keV band 
from XMM-Newton, Chandra, or Swift XRT and the 14—195 keV 
band from Swift BAT (details in Ricci et al., in prep). 


The Swift BAT 70-month all-sky survey enables us to analyze 
the optical and X-ray properties of a large sample of nearby AGN. 
It has detected 823 AGN to a flux limit of ^ 10^*^ erg cm^^ s^* 
over the whole sky ( [Baumgartner et al.|2013[ . Due to its detection 
method in the ultra hard X-rays (14—195 keV), the Swift BAT 
survey is not affected by obscuration (below Compton-thick 
levels) nor is it contaminated by stellar emission. Therefore, it also 
includes ’buried’ obscured AGN that may not have been selected 
by optical surveys (e.g., [Ueda et al.|2007). 


In this work, we use a sample of BAT AGN based on the 
70-month catalog with an optical spectroscopic sample that 
is significantly larger than previous studies (340 objects). Our 
spectroscopic analysis also includes higher ionization lines not 
used previously (i.e., [NeV]A3426) to study correlations be¬ 
tween optical line emission and X-ray emission. The sample is 
presented in Section 2 and our analysis methods are described 
in Section 3. We study the correlation between [Olll] and X-ray 
emission in Section 4.1, X-ray and other high ionization lines (i.e., 
[Ne III] A3869, He II 14686, and [Ne V] 13426) in Section 4.2, and 
the correlation between only optical emission lines in Section 4.3. 
Finally, a review of our results and a discussion is presented in 
Section 5. 


We use the AGN emission-line diagnostics (e.g., [Baldwin 


[Veilleux & Osterbrock[[1987]l as updated by [Ke 


etaLjfIMT 

ley et al. ~([2006 1 to further classify our sample. We classify as 


low ionization nuclear emission-line region (LINER) any object 
that shows a LINER classification in either the [O I] 16300 or 
[S II] 16717,6731 diagnostic. For the remaining AGN, we use 
the [Nil] 16583 diagnostic, which is most sensitive to AGN, to 
differentiate between Seyfert AGN and Hll/composite/ambiguous 
(17/340). We use the presence of a broad line in Ha to differentiate 
between type 1 (217/340) and type 2 objects (116/340). A more de¬ 
tailed description of the reduction procedures of the optical spectra, 
sample selection, and overall properties of the optical spectroscopic 
sample can be found in the BASS-I paper (Koss et al., in prep). 


3 ANALYSIS OF OPTICAL SPECTRA 

We fit our sample of optical spectra using an extensive spectro¬ 
scopic analysis toolkit for astronomy, PySpecKit, which uses a 
Levenberg-Marquardt algorithm for fitting ( [Ginsburg & Mirocha] 
[201 1[ |. We implement separate methods for fitting sources with 
narrow lines only and sources with broad lines. We fit the Hj3 
(4650-5050 A), Ha (6250-6770 A), and [On] (3300-4000 A) 
spectral regions separately. For narrow-line sources, we first fit the 
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spectrum with stellar templates using the pPXF software (Penalized 
PiXel Fitting; [Cappellari & Emsellem|2004| l to remove the galaxy 
continuum and stellar absorption features that can effect line 
measurements. All emission line fits were inspected by eye to 
verify proper fitting. 

We adopt a power-law fit (1st order) to model the continuum 
and Gaussian components to model the emission lines for different 
spectral regions. For the Fla spectral region, we fit the [O I] A6300, 
[Nii]A6548, Ha, [Nii]A6583, and [Sii] A6716, and [Sii] A6731 
lines. The narrow-line widths are tied together with an allowed 
variation of ± 8 A and the systemic redshift is determined from the 
Ha line. The relative strengths of [Nll]A6548 and [Nll]A6583 
lines are fixed (at 1:2.94). For the Hj3 spectral region, we fit 
the HeIlA4686, Hj3, [OlIl]A4959, and [OlIl]A5007 lines. 
The narrow-line widths are also tied together with an allowed 
variation of ± 8 A and their central wavelengths are defined by the 
redshift of [O III] A5007. The intensity of [O III] A4959 relative to 
[O III] A5007 is fixed (at 1:2.86). For the [O II] spectral region, we 
fit the [Nev]A3346, [Nev]A3426, [Oii]A3727, [Ne III] A3869, 
and [Ne III] A 3968 lines. 

Two Gaussian components are allowed for the Ha and Hj3 
emission lines: a narrow component (FWHM < 1000 kms^*) and 
a broad component (FWHM > 1000 kms^'). The broad HjS com¬ 
ponent is only allowed if a broad Ha component is detected. In the 
case of a non-detection of the narrow Ha or [SII] A6717,6731 line 
(i.e., due to a very strong broad Ha component, a noisy spectmm 
or bad sky subtraction), we use the FWHM of [OlIl]A5007 to 
constrain the widths of the narrow lines in the Ha region. 

To estimate the continuum for the Ha complex, we use the 
wavelength regions 5800-6250 A and 6750-7000 A. For the Hj3 
complex, we use the 4660-4750 A (except around Hell) and 
5040-5200 A regions. The continuum around the [Oil] spectral 
region is usually more complicated to fit due to a non-linear shape 
or because it lies at the blue limit of the wavelength coverage. To 
fit this blue continuum, we use the region between 3300 A and 
4000 A except where the emission lines are located (± 15 A). To 
estimate errors in line fluxes and line widths, we use a Monte Carlo 
simulation that adds noise based on the error spectrum and reruns 
the fitting procedure 10 times. The flux uncertainty for the [Olll] 
emission line is typically less than 1%. To correct our line ratios 
for dust extinction, we use the narrow Balmer line ratio (Ha/H/3) 
assuming an intrinsic ratio of R=3.1 (e.g., [Ferland & O sterbrock] 
|1986| > and the |Cardelli et al.| jl989^ reddening curve when both 
the narrow Ha and HjS are detected. Figure shows an example 
spectrum and best model fit. 

4 RESULTS 

4.1 X-ray and [O III] A5007 Flux 

We first compare the X-ray flux (2—10 keV and 14—195 keV) to 
the [O III] A5007 emission line (Figure]^. The relation between 
the [O III] corrected flux and the BAT X-ray flux shows a scatter of 
(1=0.62 dex and a Pearson correlation coefficient of Rpear = 0.35. 
The result is very similar when comparing [O III] to the 2—10 keV 
flux (ct= 0.63, Rpear = 0.38). The Pearson p-value representing the 
probability of the hypothesis to have no correlation at all between 


the two X-ray fluxes and the [Olll] flux is < 10^^® because of 
the large sample size. The correlation in luminosity space is better 
(Rpgjjr*®^ = 0.64, Rpea*** = 0.63) than in flux space, as expected 
because the same multiplicative factor (47tdf) is applied to both 
the X-ray flux and the [Olll] flux to obtain the luminositieQ 
We finally compare the [Olll] A5007 flux both with (Fjq'jjj) and 
without (Fjoin]) extinction correction to the two X-ray band fluxes. 
We find no significant improvement in the scatter of the different 
relations (~0.01 dex) using the Balmer decrement corrected [O III] 
fluxes. 

Next, we compare a one-to-one relation and the linear 
regression in logarithmic space to fit our data. For both slope fits, 
we remove the four largest outliers. For the line fitting, we use the 
python routine scipy. optimize. curve_fit which provides a 
least squares approach (OLS(Y| X)) with the Levenberg-Marquardt 
gradient method for convergence. The black solid lines in the 
panels of Figure show the best one-to-one relation (which 
would imply that TJo^n] “ ^14-195)- The best fit for the linear 
regression, shown with the red lines, is Tjom] 
and oc To test if the linear regression is 

a significantly better fit, we run a E-test. The E-test compares 
the sum of the residuals squared of two different models. We 
obtain a E-value of 7.66 for the 14—195 keV flux and of 7.81 
for the 2—10 keV flux. The resulting p-values are lower than 
1% suggesting that the linear regression between the [Olll] 
luminosity and the two X-ray bands is a significantly better model 
than a one-to-one relation. However, the improvement in the 
scatter is very small (~ 0.01 dex). Therefore, we decided to use the 
one-to-one relation (the simplest model) for the rest of the analysis. 

We also use the Astronomy Survival Analysis (ASURV) 
package from the STATCODES suite of utilities developed by 
Eric Eeigelson (Feigelson 1985). We determine the correlation 
between [Olll] and X-ray luminosity in the presence of the 13 
upper limits in [Olll]. Specifically, we use the EM algorithm, 
though the Buckley James yields similar results. We find the best 
fit correlations is which is consistent 

with the least squares approach suggesting the small number 
of non-detections do not significantly change the slope of the 
correlation. 

Finally, we use the OLS bisector method of line fitting which 
is more reliable at recovering the intrinsic slope in the presence of 
uncertainties in both the X and Y variables (e.g. Isobe 1990). Table 
[T] lists the best fit linear regressions for the luminosity-luminosity 
plots using the OLS bisector method as well as the standard 
deviation for the complete sample and various sub-samples (Type 
1-1.9, Type 2, Seyferts, LINERs, and Hll/composite/ambiguous). 
The LINERS and Hll/composite objects (from the BPT clas¬ 
sification; [Baldwin et al.|[T981| jVeilleux & Osterbrockj [T987l > 
have on average lower X-ray and [Olll] luminosities than 
Seyfert galaxies: (log = 40.5; (log qZU = 40.9; 


^ di varies here from 43 Mpc to 2.2 Gpc 
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4.1.1 Correlations with obscuration and luminosity 

The two upper panels of Figure show the ratio as a 

function of X-ray luminosity. The red dots show the average trend 
and the error bars the average 1-CT scatter. First, we notice that there 
is no significant difference between the 2—10 keV (right panels) 
and the BAT (left panels) X-ray luminosity relations. There is a 
slight decrease of the [O IIl]/X-ray flux ratio when we go to higher 
X-ray luminosities, which is consistent though with being constant 
at a 3-(T level. Figure|^illustrates the luminosity dependence of the 
flux ratio using histograms. The scatter slightly decreases towards 
high X-ray luminosity. This trend is not statistically significant, and 
is probably driven by the increase of the unobscured AGN fraction 
with high luminosities (receding torus model, e.g., |Lawrence|199H 
|Simpson||1998| |2005[ |Brightman et al.|r2015| l since unobscured 
AGN exhibit lower scatter than the obscured AGN (see central 
panels of Figuresee below). The two lower panels of Figure]^ 
show a strong increase of the flux ratio as a function of the [O III] 
luminosity. 

Dust and obscuration might still play a significant role in the 
relation between optical and X-ray emission even if the optical 
emission lines are corrected for host galaxy extinction using the 
B aimer decrement. The obscuring torus scatters a portion of the 
ionizing light from the nucleus and therefore also has an indirect 
effect on the optical emission line intensities. The two central 
panels of Figure show the L|qj|jj/Lx ratio as a function of the 
integrated column density Ajj. Our results show that the trend 
of the luminosity ratio is consistent with being constant from 
very low to very high levels of obscuration. This is true for the 
intrinsic 2—10 keV emission as well as for the 14—195 keV band. 
We obtain average ratios of log T'q||jj/Li 4 _i 95 = —2.42 and 
log = —2.01. Flowever, the scatter strongly increases 

for Ah> 10^^ cm^^. We observe a scatter of CT = 0.47 for low 
obscuration AGN (Ah< 10^^ cm^^) and of a = 0.72 for high 
obscuration AGN (Ah> 10^^ cm^^) in the case of the BAT X-ray 
band. The result is similar for the 2—10 keV emission. In order 
to quantify this change in scatter, we run a Levene’s test. This 
test checks whether the two samples are drawn from distributions 
having a similar variance. We obtain a p-value = 6.52 x 10^® 
for the BAT band and of 2.35 x 10^“^ for the 2—10 keV band, 
implying a highly significant increase in the scatter from low to 
high obscuration in both X-ray bands. 


4.1.2 Slit Size and Dilution 

Due to different slit widths (0.75-3.0") used to observe the BAT 
AGN and the various redshifts of these AGN (0.01 < z <0.4), our 
spectroscopic data cover regions between 0.2 kpc to the size of the 
entire host galaxy. In addition to that, the NLR size may vary from 
one galaxy to another and the variation is strongly dependent on 
the power of the AGN (e.g., |Mor et al.|20()^|Hainline et al.|2013| >. 
This change in the size of the observed region and in the size of 
the NLR has various implications for our measurements. First, 
the region we probe with the spectroscopic observations might 
not cover the entire ionized [O III] region for luminous sources at 
low redshift. We might therefore underestimate the [O III] flux for 
these objects. Alternatively in low luminosity AGN, the region we 
probe with the slit could be much larger than the size of the NLR 


with a significant contaminating component from the host galaxy. 

We study /Lx as a function of the physical slit size in 
kpc (lower panels of Figure to study these issues. Our results 
show that the scatter slightly decreases towards larger slit sizes 
while the ratio is constant. To further quantify this trend, 

we split our sample in two based on slit sizes > 3 kpc and slit sizes 
< 3 kpc. We observe that the scatter decreases by respectively 0.16 
dex and 0.20 dex for the 14—195 keV and 2—10 keV band. How¬ 
ever, this decrease in scatter is probably driven by the increase of 
the unobscured AGN fraction with high redshift rather than simply 
larger physical slit sizes. 

4.1.3 Sample inhomogeneity 

The BASS survey gathered optical spectra from various instru¬ 
ments with different data qualities (i.e., resolution, sensitivity, 
wavelength range). This lack of homogeneity in the data could ar¬ 
tificially increase the scatter in the [O III] - X-ray relation. For this 
reason, we also analyze the SDSS sub-sample (110 AGN) because 
of their high quality in flux calibration (Figure]^. We find that the 
scatter decreases by 0.09 dex for the 14—195 keV and 0.08 dex 
for the 2—10 keV band. In order to have a quantitative result for 
this change in scatter, we run again a Levene’s test. The resulting 
p-values for the SDSS sub-sample are respectively 0.062 and 0.35 
for the BAT and 2—10 keV X-ray bands which is not statistically 
significant. The fit of the SDSS sub-sample (red line) is similar to 
the fit of the total sample (black line) for both X-ray bands. This 
implies that we do not have any statistical offset due to systematics 
(like flux calibration issue). We conclude that there is no significant 
difference between the total sample and the SDSS sub-sample. 

4.2 X-ray correlation with high and low ionization lines 

Lines produced from high ionization (>40 eV) have the advantage 
of being less contaminated by host galaxy emission than lines 
produced from low ionization (<20 eV). Moreover, these lines are 
emitted in less extended regions, which implies that we cover a 
larger portion of their emitting regions in our observations than for 
low ionization lines. Finally, we include the Ha and Hj3 emission 
lines as these lines are either not or less affected with respect to 
other lines by variations in metal abundances, ionization levels, 
and collisional effects. 

Figure illustrates the various emission line fluxes as a 
function of the two X-ray fluxes. Our data show that the correlation 
is slightly stronger with [Nelll], Hell or [Ne V] than with [Olll], 
Ha or [Sll], especially for the BAT X-ray band. Their scatter 
is smaller by about ~0.09 dex and the correlation coefficient 
increases from Rpear ~ 0.36 to Rpear ~ 0.44. However, the [Ne V] 
emission line is only detected in optically selected AGNs with very 
bright emission lines whereas the [O III] line is detected in many 
faint AGN such as the LINERs and composite objects which has 
the effect to increase the scatter. To overcome these possible biases, 
we need to compare the same samples together. This is what we 
show in Table In this table, the reported standard deviations 
and Pearson coefficients are calculated based on a sub-sample of 
objects having both the detection of the considered line and an 
[Olll] detection. In order to quantitatively investigate if the flux 
correlation of X-rays with [S II], Ha, [Ne III], He II or [Ne V] differ 
significantly from the one with [O III], we run a z-test based on the 
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two Pearson correlation coefficients (Fisher r-to-z transformation). 
The p-value is also reported in Table|^ The result of this statistical 
analysis shows that we have p-value > 0.5 for all emission lines. 
The small differences between the linear relations are therefore not 
statistically significant. 


4.3 Relationships between optical emission lines 

Figure shows the relation between [O III] and other optical 
emission lines. We see from the scatter (~0.3 dex) that these 
emission lines are tightly correlated together. The mean flux ratios, 
standard deviations and Pearson coefficients are summarized in 
Table[^ One can especially notice the [O III] - [Ne III] relation that 
exhibits a scatter of CT = 0.24 dex. The two transitions originate 
from ions with comparable ionization energies (35.1 and 41 eV, 
respectively), and from gas with comparable 
(5 X 10^ cm^^ and 7 x 10® cm^^, respectively; 

|2007|> which may explain the small scatter. 


critical densities 
Kogure & Leung 


5 SUMMARY AND DISCUSSION 

We have presented a study of emission line properties of the largest 
hard X-ray selected sample of AGN in the local Universe. Our con¬ 
clusions are: 

(i) We find very weak correlations with a large scatter (0.62- 
0.63 dex) between the two X-ray fluxes (2—10 keV and 14—195 
keV band) and the observed and dust corrected [Olll]. The ra¬ 
tios between the X-ray fluxes and the [Olll] corrected flux are 
Fi4-195 = (267 ±22) X and ±2-10 = (102±8) x 

(ii) The correlation of emission lines with high ionization levels 

([Ne III], He II, [Ne V]) with the X-ray flux is not significantly better 
than for [O III]. This is also true for the lower ionization line [S II] 
and the Ha recombination line as compared to [Olll]. While ob¬ 
jects which have a [Ne V], He II or [Ne III] detection have a smaller 
scatter (0 'jq®)'jj = 0.49, = 0.47, = 0.45), this is only 

because they have brighter emission lines that exclude the popula¬ 
tion of LINERS or X-ray detected AGN classified as Hll/comp. In 
general, all lines show a statistically similar level of scatter with the 
X-rays independent of ionization energy. 

(iii) All the optical emission lines are strongly correlated with 
[Olll]. These correlations (cT ~ 0.3) are significantly tighter than 
what we find between [O III] and the X-rays. 

(iv) The [O IIl]/X-ray flux ratio as a function of the column den¬ 
sity (Ah) is consistent with being constant. However, we see a 
strong increase in the scatter above Ah> 10^^ cm^^ (ct<22 = 0.47, 
ct> 22 = 0.72). Our results show that the scatter slightly decreases 
towards larger physical slit sizes in kpc while the ratio LJqjjjj/Lx is 
constant for the AGN in our sample (z > 0.01). 

We find that our study agrees well with past studies of X-ray 
selected AGN, but shows a weaker correlation than [O III] selected 
samples which exclude X-ray selected AGN with weak emission 
lines (e.g., LINERs and optically elusive AGN classified as 
Hll/comp). The correlation we find between the [Olll] corrected 
luminosity and BAT X-ray luminosity agrees with the 9-month 
sample of BAT AGN (e.g., Rpear = 0.41 vs. Rpe^^ ~ 0.35-0.40; 


[Winter et al.||20 10l). Our correlation result disagrees significantly 
with the 9-month sample of [Melendez et alT| ( |2008^ which com¬ 
pared 2—10 keV X-ray fluxes from the ASCA database which 
are not calculated using detailed X-ray modeling. [Panessa et al.[ 
( [2006^ calculated the correlation between the [Olll] corrected 
flux and the 2—10 keV flux from a sample of 47 nearby Seyfert 
galaxies from the Palomar optical spectroscopic survey ( [Ho et al!] 
[1995[ (. They obtained a strong Spearman correlation coefficient of 
rj = 0.78 that is also not in accordance with this work (rs = 0.39) 
most likely because their sample was optically selected and did not 
include AGN that were X-ray weak. Our results show a slightly 
higher scatter than the [O III] - X-ray relation in [Heckman et al.[ 
l [2005t (<7 = 0.51; this work: o = 0.62) whose data come from the 
RXTE slew survey (3-20 keV). 

We find no significant evidence of a weakening or strength¬ 
ening of line emission at high AGN luminosities. This is because 
of the very weak correlation and large scatter which make a linear 
model inappropriate (e.g. Hogg et al. 2010). The improvement in 
the scatter of a linear fit is very small (~ 0.01 dex) over a simple 
one-to-one relation between [Olll] and X-ray luminosity. Using 
the OLS bisector method, the best slope between the intrinsic 
[Olll] and the 14—195 keV luminosity is 1.29 ±0.05 and the 
intrinsic [O III] and 2—10 keV luminosity is 1.23 ±0.05. This is in 
accordance with previous BAT studies (e.g., 1.16 ± 0.20; [Winter[ 
[et al.[[^10[ (. Alternatively, the EM method including censured 
data finds a slope of 0.83 ± 0.06 which is consistent with past 
studies using the same line fitting method as found in [Melendez[ 
|eraL]|2008| (0.8 ±0.1) or [Panessa et~aL]|2006| (0.82 ± 0.04). In 
summary, our best fit slope is consistent with past studies, though 
the large intrinsic scatter of the relation yields very different results 
depending on the line fitting method used. 

Our analysis shows that the relation between optical emission 
lines and hard X-rays exist, but is very weak (at best C7 ~0.5 dex), 
no matter which optical emission line ([Sll], Ha, [Olll], [Nelll], 
He II, [Ne V]) or which X-ray flux (2—10 keV intrinsic or 14—195 
keV). Additionally, calibration, obscuration, and physical slit size 
play only a small role in the scatter. This scatter is much larger 
than the scatter of X-ray compared to SED fits of unobscured 
BAT AGN using the UV (e.g., 14—195 keV scatter of 0.3 dex, 
[Vasudevan et al.[|2009[ l. While our analysis focused on the hard 
X-rays (>2 keV), the soft X-rays (<2 keV) have been shown to 
be dominated by emission lines in several Seyfert AGN whose 
extended spatial profiles show morphological correlations with the 
[Olll] (e.g. Bianchi 2006, Wang et al. 2011, Koss et al. 2015). 
Thus, the soft X-rays coming from the emission line regions may 
shower a tighter correlation that is worth investigating in future 
studies. 


These results therefore strongly suggest that the majority 
of the intrinsic scatter has physical origins related to the NLR 
gas. There are several mechanisms such as covering factor of the 
narrow-line region gas (e.g., [Netzer & Laor|1993[ Noguchi et al.[ 
[2010[ (, density dependence of [O III] ([Crenshaw & Kraemer[2005' 


differing ionization parameter ([Baskin & Laor|2005|, contamina- 
tion of the narrow emission lines from star formation (e.g., |Wd| 
[et al.[[20TT| l, errors in the Balmer decrement because of different 
distribution of Balmer lines and forbidden lines, and finally SED 
shape changes with luminosity jNetzer et al.|2006T l. Another possi¬ 
bility is AGN variability, where the X-ray emission varies on very 
short time scales (sub-parsec scale) and represents the instanta- 
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neous accretion rate while the NLR which is orders of magnitude 
larger than the corona (100 parsec to kiloparsec scale), varies on 
much longer time scales and traces the average accretion rate (e.g., 
|Mushotzky et al.|1993[|Schawinski et al.|2015^ . Future studies us¬ 
ing high signal-to-noise Integral Field Unit (IFU) observations with 
wide wavelength coverage of large samples of AGN such as the 
Multi Unit Spectroscopic Explorer (MUSE; [Bacon et al.|2010^ are 
critical to map and study the importance of these different parame¬ 
ters. Einally, the Nuclear Spectroscopic Telescope Array (NuSTAR; 
|Flarrison et al.|[20T3) hard X-ray telescope is currently observing 
many BAT AGN as part of the legacy survey (Balokovic et al., in 
prep.) with a 100 times increase in sensitivity which improves esti¬ 
mates of X-ray luminosity compared to BAT. 
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Table 1. [OIII] - X-ray luminosity relations for different sub-samples 


Y 

X 

Sample 

N 

a 

b 

(7 

<X> 

<Y> 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 



all 

327 

L29±0.05 

-15±2 

0.61 

43.9 

41.4 



Type 1-1.9 

215 

1.17±0.06 

-9.6±2.6 

0.53 

44.0 

41.6 



Type 2.0 

112 

L40±0.10 

-20±4 

0.69 

43.7 

41.1 


log Li4_195 

Seyferts 

241 

1.16±0.05 

-9.3±2.4 

0.53 

43.9 

41.5 



LINERs 

17 

1.8±0.3 

-39±13 

0.70 

43.6 

40.5 



Hll/comp/ambiguous 

17 

L2±0.2 

-11±8 

0.74 

43.7 

40.9 



all 

321 

1.23±0.05 

-12±2 

0.61 

43.5 

41.4 



Type 1-1.9 

209 

LlOiO.06 

-6.5±2.5 

0.56 

43.6 

41.6 



Type 2.0 

112 

L30±0.10 

-15±4 

0.69 

43.3 

41.1 


log if 10 

Seyferts 

238 

L13±0.05 

-7.5±2.3 

0.54 

43.4 

41.5 



LINERs 

17 

L8±0.3 

-37±12 

0.62 

43.1 

40.5 



Hll/comp/ambiguous 

16 

1.0±0.2 

-4±8 

0.62 

43.2 

40.8 


Note. — (1) variable Y ; (2) variable X ; (3) AGN type : type 1-1.9 are all the objects with broad Ha detection, type 2 
are all the objects without broad detection, LINERs are selected by the [S II] and [O I] diagrams, Hll/comp/ambiguous and 
Seyferts are selected by the [N II] diagram (except for previously selected LINERs); (4) size of sample ; (5) slope with its 
1-(J uncertainty ; (6) regression intercept with its l-fj uncertainty. The equation of the linear regression is Y = aX -I- b ; (7) 
standard deviation ; (8) mean value of X ; (9) mean value of Y. 


Table 2. Relations of various optical emission lines intrinsic fluxes with the 14—195 keV emission and the comparison with 
the [0 III] - 14—195 keV relation. 


Line 

A[eV] 

N 


^Line 

Ppcilr. lOlIl] 

^^Pear. Line 

p-value 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

[0 i]A6300 


206 

0.60 dex 

0.53 dex 

0.41 

0.46 

0.54 

[Sll]A6717,6731 

10.4 

275 

0.61 dex 

0.53 dex 

0.37 

0.38 

0.89 

Ha 

13.6 

296 

0.62 dex 

0.57 dex 

0.36 

0.34 

0.78 

[Nll]A6583 

14.5 

277 

0.62 dex 

0.56 dex 

0.37 

0.36 

0.89 

[Nelll]A3968 

41.0 

144 

0.47 dex 

0.48 dex 

0.45 

0.44 

0.92 

He II A4686 

54.4 

106 

0.45 dex 

0.49 dex 

0.49 

0.42 

0.53 

[NeV]A3426 

97.1 

51 

0.49 dex 

0.49 dex 

0.42 

0.45 

0.86 


Note. — (1) optical emission line ; (2) ionization level; (3) size of the common sample in which both the Line and [O III] are 
detected ; (4) standard deviation and (5) Pearson R coefficient of the log - log T 14-195 relation ; (6) standard deviation 

and (7) Pearson R coefficient of the log - log F 14-195 relation ; ( 8 ) p-value of the null hypothesis that the two correlation 
coefficients obtained from independent parent samples are equal. 


^ scatter of the [OlII] - 14195 keV relation only for objects with [OlII] and the other emission line (1) detection to avoid 
selection biases. 

^ Pearson correlation coefficient for the same sub-sample used for (T|oin] 
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Table 3. Correlation between [O III] A5007 and other optical emission line intrinsic fluxes. 


Line 

A[eV] 

N 

<r> 

a 

Rpear 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

[0 i]A6300 


206 

0.050 

0.36 dex 

0.84 

[Sll]A6717,6731 

10.4 

275 

0.19 

0.37 dex 

0.83 

Ha 

13.6 

296 

0.438 

0.30 dex 

0.90 

[Nll]A6583 

14.5 

277 

0.30 

0.40 dex 

0.81 

[Nelll]A3968 

41.0 

138 

0.17 

0.24 dex 

0.89 

He II A4686 

54.4 

106 

0.031 

0.25 dex 

0.88 

[NeV]A3426 

97.1 

43 

0.18 

0.35 dex 

0.79 


Note. — (1) optical emission line ; (2) ionization level; (3) size of the sample ; (4) mean flux ratio ; (5) standard deviation ; 
(6) Pearson R coefficient. 
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Figure 1. Example of the fitting procedure for MCG +06-16-28 observed with the Palomar Double Spectrograph. Top left panel: Spectrum of the object with 
all the fitted emission lines. Rest frame continuum subtracted fits are shown for the [O II] (top right), H/3 (bottom left), and Ha (bottom right) spectral regions. 
The colors denote regions used for continuum (green), different Gaussian components of the fit (blue), and model fits (red) with residuals shown in the bottom 
panel. 
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Figure 2. Plots of the relation of the [OlII] A5007 emission line strength with the BAT X-ray emission (left) and with the 2—10 keV emission (right). The 
to III] — X-ray relation exhibits a large scatter in all panels (;=e0.6 dex). The black solid lines show a fit constrained hy Y o^X with the dashed lines showing 
the ±1 dex. The red lines illustrate the linear regression. For both line fits we removed the four largest outliers. The color code comes from the the [NII] BPT 
classification. White symbols do not have any BPT classification because of undetected emission lines. The AGN showing broad Ha emission are classified 
as type 1. On the top and central left panels the grey dashed lines illustrate the 14—195 keV flux limits for 50% (left line) and 90% of the sky (right line). The 
upper limits were ignored in the analysis. Top panels: observed 10 III] flux plotted against the X-ray fluxes. Central panels: extinction corrected [OlII] flux 
plotted against the X-ray fluxes. Lower panels: extinction coiTected [O III] luminosity vs. the two X-ray luminosities. 
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Figure 3. Plots of the ratio as a function of X-ray luminosity {top panels), column density {central panels), and physical slit size {lower panels). 

The color codes represent the [NII] BPT classification. White symbols do not have any BPT classification because of undetected emission lines. The red dots 
represent average binned values and the red bars show 1-sigma deviations in the bins. The first bin of the central panels represents a upper limit where the 
column density is consistent with the Galactic neutral hydrogen. The horizontal black lines on the lower panels illustrate the 3 kpc value. The scatter increases 
for Afi > 10^^ cm“^ with no change in average For slit size, the scatter decreases with increasing slit size because of a higher fraction of unobscured 

AGN at higher redshifts, but the average is constant. The upper limits were ignored in the analysis. 
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Figure 4. Histograms of the ratio binned by X-ray luminosity (top panels) and [OlII] luminosity (lower panels). The two lower panels show a 

strong positive correlation of the ratio with the [O III] luminosity. 
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log Fi 4 _i 95 keV (o^g -Cm ^ -S ^ ) log ^eV (^^g -Cm ^ -S ^ ) 


Figure 5. Plots of the [O III] A5007 extinction corrected flux vs the two X-ray fluxes. The red triangles represent the sources observed by SDSS (110 objects). 
The black line is a one-to-one fit of the total sample and the red line a one-to-one fit of the SDSS sub-sample. The scatter and the fit of the sub-sample are 
comparable to the ones of the main sample. 
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Figure 6. Plots of the relation between five optical narrow emission lines ([SlI]A6717,6731, Ha, HeIlA4686, [Nelll] A3869 and [NeV]A3426) with the 
BAT X-ray flux {left) and with the intrinsic 2—10 keV flux (right). The plots show subsamples where the particular emission line is detected. The flux units 
are in erg ■ cm~^ • 5 “^. The subplots are vertically ordered by ionization level of the emission line. The black lines are one-to-one fits of the data. The color 
code represents the [NII] BPT classification. White symbols do not have any BPT classification because of undetected emission lines. 
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Figure 7. Plots of the relation between optical narrow emission lines and the [O III] A5007 line in flux space. The plots show subsamples where the particular 
emission line is detected. The lines are all strongly correlated with [O III]. The black lines are one-to-one fits of the data. The color code represents the [N II] 
BPT classification. White symbols do not have any BPT classification because of undetected emission lines. 










